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The differential energy spectrum of the cosmic radiation from solar modulation energies up to 5 X 10^^ eV 
is correctly predicted by a recent theory of the knee and ankle which uses only one normalization point. This 
remarkable quantitative result, spanning over many decades in energy and intensity, along with the existence of 
the second knee at 6 x lO^'^ eV, is obtained assuming constant spectral indices of individual ions at the cosmic- 
ray sources and no other critical hypotheses. In this study the chemical composition of the cosmic radiation is 
evaluated as a direct consequence of the theory. The computed mean logarithmic mass exhibits a rising trend from 
1.8 to 3.0 in the range 10^^-10^^ eV, a maximum value of 3.2 at 3 x 10^^ eV, and a characteristic lightening above 
3 X 10^^ eV up to 4 X 10^* eV. All of these distinctive features are in accord with the data of many experiments. 
Two additional consequences intrinsic to the theory are qualitatively discussed: (1) some limitative bounds on 
the mechanism accelerating cosmic rays; (2) the degree of isotropy implied by the residence time of the cosmic 
rays in the Galaxy. 



1. Introduction 

A solution to the knee and ankle problem has 
been recently proposed [1-4] with a major charac- 
teristic feature consisting of a fundamental con- 
nection between the knee and the ankle mainly 
based on some data on the Milky Way (size, mag- 
netic field, gas density) and nuclear interaction 
lengths of cosmic ions in the interstellar medium. 
The knees and ankles of individual ions are gen- 
erated by galactic cosmic rays, though a tiny ex- 
tragalactic component might be admitted. The 
upper energy boundary of 5 x 10^^ eV marks the 
iron ankle [2] in the energy region of the quasi 
rectilinear propagation of iron nuclei through the 
Galaxy. 

The present study assumes that this solution 
does exist as reported elsewhere [1-4] and it 
explores some consequences regarding: (1) the 
chemical composition of the cosmic radiation 
above the knee up to 5 x 10^^ eV; (2) some basic 
constraints on the acceleration mechanism in the 
Milky Way; (3) the residence time of galactic cos- 
mic rays and its relevance pertaining to the level 
of isotropy in the arrival directions of the cosmic 
rays at Earth. 

No new elements of the theory are introduced 



in this study. 

Other salient consequences of the common ori- 
gin of the knee and the ankle are not discussed 
due to the paper size limitation. For example, the 
existence of the second knee around 6 x 10^^ eV 
[5,6] is an intrinsic property of the energy spec- 
tra resulting from this theory. A second example 
is that the iron knee is less pronounced than the 
knees of the lighter ions, in particular those of 
helium and proton. This circumstance may be 
quite fruitful while searching for the iron knee 
in the energy region 

10i7_ioi8 eV as aimed at 
by the Kascade-Grande experiment, presently in 
data taking [7] . 

The quantitative agreement of this theory with 
the experimental data of Akeno [5] , Haverah Park 
[8] and Auger [9] is illustrated in figure 1 (other 
data are in fig. 11 and 12 of ref.2). Note that in 
figure 1 only one energy point at 10^* eV is used 
to normalize the predicted cosmic-ray intensity. 

The energy spectra of individual ions or group 
of elements are also in agreement with the exper- 
imental data as reported elsewhere (ref.4, fig. 5). 

Notice that the agreement of the computed 
spectrum with data displayed in figure 1 does not 
critically depend on the ion abundances nor on 
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Figure 1. The thick hne is the result of calcu- 
lation with the differential cosmic ray intensity 
normalized to l.f 6 x 10"^* particles/ (m^ sr s eV) 
at 10^"' eV. The data are from the Akeno [5], Hav- 
erah Park [8] and Auger [9] experiments. 

the specific normalization energy 0. The accord 
persists, whatever the input parameters, provided 
that they are comprised within the observational 
limits defining the Galaxy and those defining nu- 
clear interaction cross sections. 

2. Basic elements of this theory related to 
the present calculation 

A logical partition of this theory would iden- 
tify two areas: (1) the simulation of cosmic-ray 
trajectories in the Galaxy, detailed in the quoted 
papers [1-4] and others [10-12]; (2) the introduc- 
tion of a new postulate, a fundamental part of the 
theory, which states that the spectral indices of 
all ions have constant, common values of about 

^Michelangelo Ambrosio and Carla Aramo from Univer- 
sita di Napoli, Italia, suggested to us to normalize the 
differential energy spectrum at any energy in the interval 
IQli-lo" eV where numerous balloon and satellite data 
are available. This normalization choice, unlike that at 
lO^S and 10^^ eV [2], besides its superior reliability and 
precision, makes more impressive the agreement between 
theory and data due to the larger excursion between the 
normalization energy (10^'^ eV) and the extreme end of 
the predicted spectrum (5 X 10^^ eV). 
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Figure 2. Energy spectra of 11 individual ions 
measured by balloon-borne and satellite experi- 
ments suggesting and founding the postulate of 
the constant spectral indices. 

2.65 up to 5 X IQis eV. 

The postulate is well anchored to numerous ob- 
servations gathered by balloon and satellite ex- 
periments in the energy interval 10^^-10^^ eV. 
Figure 2 shows the energy spectra of individual 
ions: proton, helium, carbon, oxygen, neon, mag- 
nesiun, silicon, sulfur, argon, calcium and iron. 
The data plainly indicate that all spectra are 
straight lines in logarithmic scales and that the 
indices accumulate around a common value of 
about 2.65 as shown in figure 3. 

Even more important than balloon data is the 
Fe spectrum measured by the Kascade Collabo- 
ration [13] in the attempt to detect a signal of 
the yet unobserved Fe knee. In the large energy 
interval 10^^-8 x 10"'^^ eV, a constant index close 
to 2.64 has been measured. Thus, the Fe spec- 
trum measured by Kascade enlarges to 6 energy 
decades the interval where the values of the spec- 
tral indices are constant and almost equal. This 
fundamental observation of the Kascade Collab- 
oration reinforces the basis of the postulate. 

The intensity of the computed cosmic radiation 
can be normalized to experimental data at any 
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Figure 3. Spectral indices of cosmic ions from 
Hydrogen to Nickel measured by balloon-borne 
and satellite experiments at 10^^ eV (filled dots) 
and 10^'' eV (open dots). All indices extracted 
from the spectra around 10^^ eV fall between 2.6 
and 2.7. 

arbitrary energy in the interval 10^"'^-5 x 10^'' cV. 

The set of parameters associated to the ion 
groups is concisely termed ion blend. The present 
study adopts two ion blends: one compiled by 
Wiebel-Sooth [14] at lO^^ eV (LE blend, for Low 
Energy) and an analogous compilation [15] made 
by us at 10"^^ cV (HE blend, for High Energy) 
which includes recent balloon experiments at high 
energy [16-17]. Table 1 reports the parameters of 
these two ion blends. The relevant differences be- 
tween HE and LE blends reside, respectively, in 
the proton index, 2.67 instead of 2.77 and a larger 
proton fraction, 0.42 instead of 0.34. A paucity of 
elements heavier than helium by about 1-2 per- 
cent characterizes the LE blend compared to the 
HE blend. 

3. The chemical composition of the cosmic 
radiation above 10^^ eV 

Injecting uniformly in the galactic disk the 
same amounts of ions at given energies in the in- 
terval 10^^-5 10^® eV, the ion intensities in the 
local galactic zone (the Solar system) take the 
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Figure 4. Relative intensities of cosmic rays 
reaching the local galactic zone computed by 
assuming a uniform spatial distribution of the 
sources in the disc volume, equal abundances 
of the ions and (artificial) flat spectral indices. 
These last two constraints simplify the analysis 
at an intermediate step; they are replaced by the 
correct empirical parameters termed ion blends, 
specifying ion abundances and spectral indices. 
These curves represent the end product of the 
simulation of the cosmic-ray trajectories in the 
Galaxy. 

forms reported in figure 4. These spectra are the 
result of the first of the theory in the very 
schematic partition mentioned in Section 2. 

Any features of these spectra are thoroughly 
determined by the theory: the difference between 
the high plateau and the low plateau (see Sec- 
tion 9 in ref.l), the hierarchy in the ion intensities 
which is controlled by nuclear cross sections and 
the energy regions where the steep descents of 
the intensities (the knees) occur. The low-energy 
sides of the energy regions where cosmic ions pro- 
ceed almost unbent in the Galaxy correspond to 
the regime of quasi rectilinear propagation and 
they mark the ankles of individual ions. 

Figure 4 shows that the ion propagation 
through the Galaxy leaves almost unchanged the 
indices defined at the cosmic-ray sources. The 
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Figure 5. Computed ion spectra for the LE ion 
blend (10^^ eV). The thick hne is the complete 
spectrum. 

variations amount to 0.06 for protons and to 0.04 
for Fe nuclei in the interval lO^'^-lO^^ eV. 

When the spectral indices and ion abundances 
are included in the calculation, the regular, 
smooth spectra of figure 4 spread in intensity as 
displayed in figures 5 and 6 for the two ion blends. 
The upper curve in figure 6 (and also in figure 5) 
is the complete spectrum resulting from the sum 
of partial spectra of individual ions or groups of 
ions. 

The features of the energy spectra of two ex- 
treme ions, proton and Fe, are similar to those of 
the intermediate ions and, accordingly, the sub- 
sequent discussion is limited to protons and Fe 
nuclei. Consider for example the partial spec- 
tra of figure 6 (HE blend). The steep descent 
of the proton spectrum below 10^^ eV (the pro- 
ton knee) caused by galactic phenomena is trans- 
formed above 10^* eV into a much harder spec- 
trum controlled only by the observed index of 2.67 
of Table 1. The Fe spectrum has the same be- 
haviour but it extends, approximately, to an en- 
ergy decade higher than that of the proton. 

The intensity gaps between the high and low 
plateaux in the proton and Fe spectra of figure 
6, clearly visible in figure 4, dominate the fea- 
tures of the chemical composition in this energy 
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Figure 6. Computed ion spectra for the HE ion 
blend (lO^** eV). The thick line is the complete 
spectrum. 

region. Because of these characteristics, there is 
an evident increase of the proton fraction rela- 
tive to the Fe fraction above 3 x lO^'^ eV. This 
trait of the energy spectra becomes evident when 
the relative amounts of light and heavy ions are 
isolated as displayed in figure 7. The light com- 
ponent (proton and He) has a minimum around 
3 X 10^'' eV, then it increases regaining its status 
of a large component of about 34 percent above 
fO^^ eV. Above 3 x 10^"^ eV a moderate lighten- 
ing of the cosmic-ray spectrum takes place with a 
chemical composition dominated by a heavy com- 
ponent. A similar behaviour is exhibited by the 
spectra of figure 5 (LE blend) . 

Above 10^^ eV the proton and He fraction of 
the LE blend (light ions in figure 7) is 30 percent 
while the Fe and Ca fraction (heavy ions in fig- 
ure 7) is about 50 percent at 10^^ eV. The light 
and the heavy fractions of the HE blend become 
almost equal at 5 x 10"'^^ eV amounting to 34 per- 
cent. These significant differences in the chem- 
ical compositions resulting from the two blends 
are due to the softer spectral index of the pro- 
ton, 2.77, for the LE blend instead of 2.67 for the 
HE blend, since the Fe indices are almost equal 
in both blends H. 



theory of particle acceleration in pulsar atmospheres 
predicting constant spectral indices up to 5 X 10^® eV was 
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Figure 7. Relative amounts of light and heavy 
ions versus energy for the LE (thin line) and HE 
(thick line) ion blends. Notice the minimum of 
the light component aroimd 3 x 10^^ cV in both 
ion blends, followed by the lightening of the chem- 
ical composition. 

4. Comparison of the computed <ln{A)> 
with the experimental data 

Presently, experiments operating above 10^^ 
eV cannot distinguish individual ions or a re- 
stricted group of elements as accomplished, for 
example, by the Kascade experiment below 0.8 x 
10^^ eV. At very high energy, experiments deter- 
mine the mean logarithmic mass of the cosmic 
radiation, <ln{A)>, a physical quantity whose 
measurement exploits the properties of the gi- 
ant atmospheric cascades requiring long elabora- 
tions. Around the knee energy region the cosmic- 
ray composition is extracted from the relative 
amounts of the muon and electron component of 
the atmospheric cascades at ground (see, for ex- 
ample, Eas-top and Macro data [18]). Above lO^'^ 
eV experiments measure the atmospheric depth of 
the shower maximum, X^ax^ which is converted 
into <ln{A)> using a procedure described else- 



conceived by Guido Pizzella in 1970 [19]. Some objections 
raised to the diffusive slioclc acceleration mechanism men- 
tioned in the next Section 5 can also be applied to this 
theory. 



< 




Ion blend at 10 :eV 
Ion blend at 10"eV 



10^^ 10^^ 10'^ 10^^ 10^^ lo'^ 10^^ 



10^^ 10^° 



Energy (eV) 

Figure 8. Mean logarithmic mass of the cosmic 
radiation resulting from this theory of the knee 
and ankle in the energy interval 10^^-5 x 10^^ eV 
for both ion blends. 

where [23a]. 

Figure 8 gives the mean logarithmic mass ver- 
sus energy of the galactic cosmic rays resulting 
from the spectra shown in figures 5 and 6. 

The highly distinctive feature of the computed 
chemical composition (i.e. <ln{A)>) above 3 x 
10^^ eV is a definite lightening which occurs af- 
ter a constant increase at lower energies, in the 
band IQi^-lOi^ eV. Above 5 x 10^^ eV the chem- 
ical composition becomes almost constant up to 
the maximum energy considered here, 5 x 10"'^^ 
eV. Note that the inversion energy at 3 x 10^^ 
eV (HE blend), where an heavy composition be- 
comes lighter, is rather insensitive to the ion 
blend (5 x lO^'^ cV for the LE blend). The char- 
acteristic energy of 3 X 10^^ eV of the maximum 
is not a trivial result of the theory since the four 
peculiarities of <ln{A)>, e.g. (1) increasing from 
the preknee value of 1.6, (2) reaching a maximum, 
(3) decreasing, (4) and finally, levelling off to the 
value of 2.6, might have occured in another re- 
gion, still preserving all of the features (l)-(4). 

The comparison of the computed <ln{A)> 
with the data of three representative experiments 
is made below in the range 10^^-5 x 10^® eV. 

The <ln{A)> from Kascade data [13] is shown 
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Figure 9. Comparison of tlic coniputcd moan log- 
arithmic mass with Kascade [13], Yakutsk [20] 
and Fly's Eye [21] data (selected out of many oth- 
ers). 

in figure 9 along with that extracted from 
Yakutsk [20] and Fly's Eye [21] data. Calcula- 
tion and data agree exhibiting a common rise of 
the mean logarithmic mass from about 1.8 up to 
3.0 in the energy band 2 x lO^^-O.S x 10^^ eV. No- 
tice that both algorithms adopted in the Kascade 
data analysis (QGSjet and Sibyll) exhibit a rising 
trend of <ln{A)> versus energy. The Kascade 
data (open triangles) elaborated by the QGSjet 
algorithm at 2-3 x 10^^ eV are inconsistent with 
the extrapolated balloon data and hence excluded 
in the subsequent figure 10. Since individual ion 
spectra measured by the Kascade experiment are 
in accord with those calculations, the agreement 
on <ln{A)> in figure 9 is not surprising. No- 
tice that the Kascade ion spectra obtained by the 
unfolding algorithm [22] , which is independent of 
the hadronic interaction models, give a <ln{A)> 
in good accord with this theory. 

The Yakutsk data cover both the knee and an- 
kle energy regions. A clear rise of <ln{A)> be- 
tween 10^^-10^^ eV is observed reaching maxi- 
mum values of about 2.75 in the band 10^^-10^^ 
eV, then above 10^^ eV, a decreasing trend dom- 
inates. Figure 9 indicates that the Yakutsk data 
are in accord with those from Kascade, in the lim- 
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Figure 10. Comparison of the computed <ln{A)> 
with the Tunka [24], Kascade [13] and Auger 
[23] data. A significant discrepancy between this 
calculation and the Auger data appears above 
4 X 10^^ eV. The data represented by crosses are 
obtained by decreasing by 5 percent the X^ax 
measured by the Auger Collaboration [26] , which 
amounts to 35 g/cm?. This arbitrary operation 
would force an agreement with the theory in the 
band 4 x x 10^^ eV. 

ited energy interval where data overlap, and with 
this theory. Fly's Eye data in figure 9 show a 
<ln{A)> thoroughly consistent with the theoret- 
ical shape reported in figure 8 in the sense that 
there is an increase followed by a maximum, and 
finally, a decreasing trend. Unlike Yakutsk data, 
which have a broad maximum over two energy 
decades, the maximum of <ln{A)> in Fly's Eye 
data occurs in a narrower energy band, 2-6 x 10^'' 
eV, well consistent with these calculations. 

A clear discrepancy exists between the com- 
puted <ln{A)> and that extracted from the 
Auger data above 4 x 10^^ eV as shown in fig- 
ure 10. A similar discrepancy also appears with 
the HiRes data, shown in the subsequent figure 
13-a, which are well below this calculation. 

In this paper all the values of <ln{A)> ex- 
tracted from the measured Xmax use the Heck 
curves shown in figure 11 [25a]. 

Figure 11 reports X^ax versus energy for pro- 
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Figure 11. Atmospheric depths of shower maxima 
of proton and Fe nuclei versus energy. The upper 
and bottom double curves have been computed 
by Heck [25a] for two different hadronic interac- 
tion models. The X^ax resulting from this theory 
of the knee and the ankle (thick line inbctwcen) 
is compared with Tunka, Yakutsk and Fly's Eye 
data. 



ton and iron with two hadronic interaction mod- 
els as calculated elsewhere by others [25a]. The 
curves of figure 11 in the range 3 x 10^^-5 x 10^® 
eV do not substantially differ from those of the 
Auger Collaboration [23]. 

Figure 12 reports the X^ax measured by the 
Auger experiment [23a] corresponding to the 
same data of <ln{A)> shown in figure 10 in or- 
der to restate the discrepancy between this cal- 
culation and the Auger data. In the low energy 
interval are also shown for comparison the data 
of Spase2/ Vulcan [27] and Casa-Blanca [28]. 

The available data on chemical composition 
above 10^^ eV in different experiments are not 
coherent and some acute inconsistencies denoted 
here a, (3, 7, 5, e and Q emerge: 

(a) the mean logarithmic mass of about 3-4 
at about 0.8 x 10^'^ eV determined by Kascade 
shown in figure 9 could not match that of 1.7-1.9 
inferred by us from the Auger data. The Kas- 
cade data ends at 0.88 x 10^^ eV while the Auger 
data initiates at 4.45 x 10^^ eV. It seems unnat- 



1000 



900 - 



800 



S 600 



500 



400 



300 



200 



- Au9er(2007) 

: ★ CASA-BLANCA 

- A SPASE2-VULCA 


N(1998) 


































Proton 


















5^ 


Fe 










Present calculation 










Balloon oata 
QGSjet (Hec 
SIBYLL (Hec 


I) 









10 



10 



10 



10 10 10 

Energie (eV) 

Figure 12. Atmospheric depths of shower max- 
ima of proton and Fe nuclei versus energy. The 
Xniax resulting from this theory of the knee and 
the ankle (thick line inbetween) is compared with 
Auger, Spase2/Vulcan and Casa-BIanca data. 



ural that an abrupt change in the chemical com- 
position would take place in such a small energy 
band 1-4 x 10^'^ eV. Note that <ln{A)> from Kas- 
cade data is extracted using individual ion spectra 
and not X^ax- Since an increase of 1.5 units in 
<ln{A)> is independent of the hadronic models 
used by the Kascade experiment and since bal- 
loon data have a <ln{A)> of about 1.75 in the 
range 10^-10^^ eV, a value of 3.2 around 10^'^ is 
a solid, initial reference for the experiments oper- 
ating above 10^^ eV. 

[fi) Auger data on <ln{A)> shown in figure 10, 
duo to the small error bars, unambiguously sepa- 
rate two energy regions where the chemical com- 
position has opposite trends with energy: below 
2 x 10^^ cV decreases with energy while above 
2 X 10^^ eV increases, reaching the mean value 
of about 3 with a large error bar (the last data 
point) . This rising trend badly contrasts from the 
opposite trend of Yakutsk and Fly's Eye data in 
the interval 10^*-4 x 10^^ eV, and also differs from 
an approximate flat trend of the HiRes data in the 
same energy band. The focus is on the trend of 
<ln{A)> versus energy since absolute values are 
difficult to be ascertained due to the inaccuracies 
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Figure 13. Comparison of the computed mean 

logarithmic mass with the Kascadc [13], Hires 
[29] and Haverah Park [30] data. There is a 
large discrepancy between this calculation and 
the Hires data. The Haverah Park data of figure 
13-b, though mimic the bell form of the computed 
curves, they substantially differ from the present 
calculation. 



of the hadronic interaction models which affect 
versus energy [26] . 

(7) Comparing Yakutsk and Fly's Eye data on 
<ln{A)> at 2 X 10^^ eV (on average, about 2.2) 
with those of the Auger experiment (1.12 ± 0.15) 
a large gap is registered. 

((5) Shifting rigidly the Auger data on <ln{A)> 
in such a way to join the corresponding theoret- 
ical values of about 2.7 around 2 x IQ-'^* eV of 
the HE blend (figure 10), the resulting <ln{A)> 
at 4.2 X 10^^ eV would be 4.65 ± 0.50, a rather 
surprising high value implying a cosmic radiation 
devoid of proton and He nuclei. If the <ln{A)> of 
2.2 around 2x10^^ eV, resulting from the Yakutsk 
and Fly's Eye experiments, is used in the shift (in- 



stead of the theoretical value of 2.7), the Auger 
data at 4.2 x 10^^ eV would still have an ultra- 
heavy chemical composition. This bland form of 
reductio ad absurdum ventilates that the observed 
Xmax versus energy of the Auger experiment not 
only would necessitate a rigid shift but also a sub- 
stantial deformation of the entire curve. 

(e) Notice that the simulated X„iax distribu- 
tion of the Auger experiment around 10^^ eV has 
a width (FWHM) of about 104 g/cw? for pro- 
tons and 44 g/crn? for Fe nuclei, while the corre- 
sponding Xrnax are, respectively, 750 gjcrr? and 
680 gjcrr? [26b]. From a model- free analysis of 
the Auger Collaboration an X^ax resolution of 20 
g/cw? is estimated [23b]. The uncertainty of 22% 
[23a] in the energy scale is equivalent to about 5 
g/cm? taking into account the slopes of Xmax 
versus energy (fig. 3, ref.23a). Finally, the differ- 
ent interaction models globally contribute to the 
uncertainty with about 15 g/cm?. 

Suppose to arbitrarily and rigidly shift by 
20 gjcrr? the X^ax versus energy observed by 
the Auger experiment. Since the correspond- 
ing <Zn(A)> resulting from this arbitrary shift 
at 2 X 10^^ eV cannot decrease downwardly (in 
such an hypothetical situation the cosmic radia- 
tion is devoid of all nuclei heavier than He), the 
<ln{^A)> must displace upwardly in figure 10. 
The <lniyA)> resulting from this artificial shift 
would imply smaller atmospheric depths, i.e. an 
overabundance of heavy nuclei, which in turn en- 
tails a width of the X^ax distribution narrower 
than that expected from a cosmic radiation satu- 
rated by protons and He nuclei. The conclusion 
is that the atmospheric cascades in the present 
Auger data sample should exhibit a mismatch be- 
tween (j{Xmax) and Xmax- 

(C) Figure 13 reports the <ln{A)> extracted 
from the Hires data [29] and that from Haverah 
Park experiment with a revised data sample [30]. 
The conversion of the measured proton fraction 
Fp (see fig. 11, ref.[30]) to <ln{A)> is obtained 
by assuming that the rest of cosmic-ray fraction, 
1 — Fp, consists of Fe ions. The <ln{A)> of the 
Hires experiment has a decreasing trend in the 
band (1 - 8) x lO" eV while the Haverah Park 
experiment has the opposite trend in the same en- 
ergy interval. The average difference in the values 



of <ln{A)> in the two experiments amounts to 
about 1.4 units of <ln{A)> in the range 5 x 10^^- 
10^^ eV, a large discrepancy. 

5. Some bounds on the unknown engine 
accelerating galactic cosmic rays 

The engine accelerating cosmic rays in the 
Milky Way will release them ubiquitously with 
constant spectral indices, in the limited energy 
range explored here 10^^-5 x 10^^ eV. This is a 
basic tenet of this theory of the knee and the an- 
kle. An outdated, still widespread hypothesis dic- 
tates that the bulk of cosmic rays are accelerated 
in supernovae remnants through the mechanism 
of diffusive shock acceleration. How does this hy- 
pothesis compare with the postulate of constant 
spectral indices ? The comparison is immediate. 

Empirically, it has been demonstrated in a set 
of 12 supernova remnants that the maximum en- 
ergies of the electrons rarely cixceed 10^^ eV [31]. 
Protons should conform to this limit. 

Theoretically, particles would attain a maxi- 
mum characteristic energy by the diffusive shocks 
acceleration in supernovae remnants due to the 
limited number of kicks they can acquire during 
the supernova remnants lifetimes. For example, 
in a classical paper [32] the maximum proton en- 
ergy is estimated to be around 10^^ eV and ex- 
tended to lO^'' eV in particular environments [33]. 

If hadronic cosmic rays were accelerated in su- 
pernovae remnants with a maximum limit in en- 
ergy of 10" eV (or 10^^ eV), it would follow that 
another engine, or more engines, would have to 
operate the acceleration in the huge energy inter- 
val, 10"'^'*-5 X 10^^ eV in order to comply with the 
postulate of constant spectral indices supported 
by the data as shown in figures 2 and 3. The 
fact that the spectrum of cosmic rays in logarith- 
mic scales of intensity and energy is an approxi- 
mate straight line for many orders of magnitudes 
woiild imply a perfect matching of different accel- 
erating engines, in different energy bands, active 
in cosmic regions of dissimilar sizes, having parti- 
cle injection efficiencies timed to each engine. But 
the matching of these heterogeneous hypothetical 
engines leading to constant indices in the range 
10^^-5 X 10^® eV appears implausible. 
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Figure 14. Radial distribution of galactic cos- 
mic rays (open dots) [36,37] producing the 7 
rays with energy greater than 100 MeV as de- 
tected by the EGRET experiment. The ra- 
dial distribution of supernovae remnants in the 
Milky Way (solid lines) adjusted to the func- 
tion r exp{—r/ro), where ro=3.7 kpc [38] normal- 
ized to the cosmic-ray intensity extracted from 
EGRET data at r = 5 and r = 15 kpc. 

In the above reasoning it is silently implied that 
the diffusive shocks would accelerate cosmic rays 
below 10^** eV. The same argument applied in the 
low energy end of the cosmic ray spectrum, be- 
low 10^^ eV, taking into account the data shown 
in figures 2 and 3, relegates the diffusive shock 
acceleration in the domain of imagination. 

This inference is not isolated because a num- 
ber of empirical objections, from different ob- 
servational areas, constrast with the hypothesis 
of diffusive shock acceleration in supernova rem- 
nants: (1) since the total power transported by 
galactic cosmic rays is 10^^ erg per second [11] 
and the average kinetic; energy converted into the 
bulk motion of the supernovae is 10^^ erg, occur- 
ring at a rate of 2-3 explosions per century, the 
resulting conversion efficiency for cosmic-ray ac- 
celeration is too high compared to that observed 
in other astrophysical environments. For exam- 
ple, in nova explosions the acceleration efficiency 



of relativistic electrons is estimated to be 1% [34]. 
(2) Direct observations of 7 ray intensities from 
many supernovae remnants (for example Cas A) 
indicate a deficit of primary protons [35]. (3) The 
radial distribution of supernova remnants is not 
compatible with the radial distribution of cosmic 
rays producing the observed gamma ray spectra 
[36,37]. Figure 12 compares the two radial distri- 
butions normalized to the data points at the ar- 
bitrary galactocentric radii of 5 and 15 kpc where 
the common values are, respectively, 1.2 and 0.50 
[38]. Note that the prominence of the inconsis- 
tency results from the relative high precision mea- 
surements of the two data samples. Any of the 
two normalizations or others (not shown) equally 
evince disagreement. 

Note that the observational evidence that elec- 
trons of 10^^ eV are detected in supernova rem- 
nants is not a label that the diffusive shock accel- 
eration does actually accelerate them. Plausibly, 
an ubiquitous mechanism, still unknown and still 
undescribed in detail, operates in supernova rem- 
nants and in other astrophysical sites as suggested 
by 7 ray spectra observed in the range lO^^-lO^'' 
eV. The 7 ray spectra, having indices in the range 
1.8-2.4 as detected by Hess, Cangaro, Veritas and 
other observatories, indicate that the processes 
generating 7 rays in supernova remnants are the 
same at work in 0-B star associations, pulsar 
atmospheres, giant molecular clouds and galac- 
tic nuclei. These findings directly indicate that 
cosmic-ray acceleration is pervasive. 

The conclusion is that the acceleration of both 
hadronic cosmic rays and electrons through the 
diffusive shock acceleration in supernovae rem- 
nants is not compatible with the observed con- 
stant spectral indices and hence unreal 



^Here it is suggested that the mechanism of diffusive shock 
acceleration does not operate in nature the high energy ac- 
celeration of the bulk of the cosmic rays. The cosmic-ray 
acceleration at high energy in the magnetic clouds wan- 
dering in the Galaxy, proposed by Enrico Fermi in 1954, 
has been proved not to operate in nature [39] and it consti- 
tutes another example of failure of historical importance 
in the quest for the true mechanism accelerating cosmic 
rays. 



6. Residence time of the cosmic rays in the 
Galaxy and anisotropy 

The explanation of the knee and ankle [1-4] im- 
plies a residence time of galactic ions with highly 
distinctive features. The residence time is directly 
proportional to the grammage because at high en- 
ergy all ions travel at the same velocity, the speed 
of light. 

In the following, using qualitative arguments, 
it is suggested that the residence time implicit 
to this theory of the knee and the ankle can be 
consistent with the isotropy level of the cosmic 
radiation measured by many experiments at low 
energies, i.e. below 10^^ eV. 

A long-standing, unsolved, severe problem in 
Cosmic Ray Physics is the large discrepancy be- 
tween the computed anisotropy level and the 
measured one. Particle propagation processes 
through the Galaxy, computed by analytical 
methods which assume diffusive propagation, de- 
termine the anisotropy in the arrival direction of 
cosmic rays at Earth. The computed anisotropy 
at 10^** eV [40] is an order of magnitude higher 
than that measured. 

A similar evaluation with a different computa- 
tional method has been accomplished by Hillas 
[41,42] adopting a residence time of the form 
T^K/E\ where K is a normalization constant, 
E the cosmic-ray energy and 5 a suitable index. 

From low energy nuclear fragmentation data 5 
was believed to be in the range 0.5-0.6, while 
adopting some forms of theoretical turbulence 
(Kolmogorov and Kraichnan) in the interstellar 
medium 5 may decrease to 0.33. The lower the 
value of 5, the longer the residence time of cos- 
mic rays in the Galaxy. Longer residence times 
necessarily imply higher level of isotropisation. 

The residence times intrinsic to this theory of 
the knee and the ankle for He and Fe ions be- 
low 10^^ eV have been calculated. In the range 
10"-10i^ eV 5 results 0.08 [43], while in the range 
IQiS-io" eV (5==0.35. Due to the smaU value of 5 
below 10^^ eV an isotropisation level higher than 
that resulting from the quoted calculation [42], 
which adopted (5=0.6 and (5=0.33, is expected. 
Hence, by this indirect and involved argument, 
the anisotropy expected from this theory below 



10^ eV should be closer to observations than that 
of the quoted calculation [42] . 

Above 10^^ eV a direct calculation has not yet 
been performed due to the large number of events 
required for the evaluation of the isotropy level. 

7. Conclusions. 

Data on <ln{A)> are in accord with the 
present calculation about: (1) the rising trend of 
<ln{A)> in the band 10^^ to 3 x 10^^ eV observed 
in many experiments. (2) The maximum absolute 
value of about 3.2 of <ln{A)> around 3 x 10^"^ eV 
at the confluence of Kascade [13], Yakutsk [20] 
and Fly's Eye [21] data. (3) The lightening of 
<ln{A)> in the interval 3 x lO" to 4 x 10^® eV. 
(4) The particular energy of 3 x 10"'^'^ eV where the 
heavy cosmic-ray composition initiates a charac- 
teristic decline. (5) The relative lightening of the 
cosmic radiation in the band 4 x 10^^ to 2 x 10^^ 
measured by the Auger experiment, which is in 
accord with this calculation, though its absolute 
value is incompatible (fig. 10). 

The Hires data shown in figure 13-a completely 
disagree with the present calculation based on 
spectral indices of Table 1. 

The fourfold circumstance 1, 2, 3 and 4 men- 
tioned above and the existence of the second knee 
(see figures 5 and 6) constitute an unmistakable 
multiple fingerprint of this theory anchored to the 
available experimental data. 

The failure of the present calculation above 
10^* eV, explicit in figures 10 and 13-b, does not 
affect the theory of the knee and the ankle [1- 
4] but only the ion indices of Table 1. The Fe 
index of 2.64 as measured by the Kascade Col- 
laboration, instead of 2.62 of Table 1, and simi- 
lar changes at less than 1% level, suffice to rec- 
oncile the theory with the experimental data on 
<ln{A)> above 10^^ eV still preserving the ac- 
cord below this energy. 

The light chemical composition resulting from 
the Auger experiment above 10^^ eV, due to 
its superior precision compared to other exper- 
iments, if confirmed and assessed, would signal 
a new fascinating aspect of the cosmic radiation, 
absolutely simple and elegant: the abundances of 
cosmic ions at the sources observed at low ener- 



gies remain unchanged up to very high energies, 
5 X 10^^ eV. The upheaval in the ion abundances 
registered by all experiments between 10^^ and 
10^* eV is caused by some galactic characteristics 
which, surprisingly, tolerate the pervasive mech- 
anism accelerating cosmic ions in the Milky Way. 
This image of the Nature is a vivid echo of the em- 
pirical evidence that the spectral index of cosmic 
radiation of 2.7 below 10^^ eV persists unaltered 
in the huge energy interval 4 x 10^^ to 5 x 10"'^^ 
eV. 

Table 1 

Parameters of this theory of the knee and ankle 
termed Low Energy (LE) and High Energy (HE) 
ion blends. Intensities are E'^'^ x flux in units 
of {m^'^sr^^s^^eV^-^). 



Blend. LE HE 

10^^ eV 10^^ eV 

% 7 % 7 

H 42.4 2.77 33.5 2.67 

He 26.5 2.64 27.3 2.64 

CNO 11.9 2.68 13.6 2.65 

Ne-S 9.2 2.67 10.9 2.63 

Ca(17-20) 1.2 2.67 2.8 2.63 

Fe(21-28) 8.7 2.59 12.0 2.62 

Flux 7 Flux 7 

H 1.15 10" 2?77 3.93 lO^" Z67 

He 7.19 10^6 2.64 3.20 lO^*' 2.64 

CNO 3.24 10^6 2.68 1.60 lO^*' 2.65 

Ne-S 2.50 10^6 2.67 1.28 lO^^ 2.63 

Ca(17-20) 3.14 lO^^ 2.67 3.27 lO^^ 2.63 

Fe(21-28) 2.36 10^^ 2.59 1.41 10^^ 2.62 

Tbtai 2.71 10" 2?70 1.18 10" 2M 
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